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We study a large-scale instability in a sheared nonhelical turbulence that causes generation of 
large-scale vorticity. Three types of the background large-scale flows are considered, i.e., the Couette 
and Poiseuille flows in a small-scale homogeneous turbulence, and the "log- linear" velocity shear 
in an inhomogeneous turbulence. It is known that laminar plane Couette flow and antisymmetric 
mode of laminar plane Poiseuille flow are stable with respect to small perturbations for any Reynolds 
numbers. We demonstrate that in a small-scale turbulence under certain conditions the large-scale 
Couette and Poiseuille flows are unstable due to the large-scale instability. This instability causes 
formation of large-scale vortical structures stretched along the mean sheared velocity. The growth 
rate of the large-scale instability for the "log- linear" velocity shear is much larger than that for 
the Couette and Poiseuille background flows. We have found a turbulent analogue of the ToUmien- 
Schlichting waves in a small-scale sheared turbulence. A mechanism of excitation of turbulent 
Tollmien-Schlichting waves is associated with a combined effect of the turbulent Reynolds stress- 
induced generation of perturbations of the mean vorticity and the background sheared motions. 
These waves can be excited even in a plane Couette flow imposed on a small-scale turbulence when 
perturbations of mean velocity depend on three spatial coordinates. The energy of these waves is 
supplied by the small-scale sheared turbulence. 

PACS numbers: 47.27.N-; 47.27.nd 



I. INTRODUCTION 

Large-scale vortical structures are universal features 
observed in geophysical, astrophysical and laboratory 
flows (see, e.g., [H, 0, [3|, y, [sl, [g])- Formation of vorti- 
cal structures is related to the Prandtl secondary flows 
(see, e.g., 0, S H, [l^). A lateral stretching (or "skew- 
ing" ) by an existing shear generates streamwise vorticity 
that results in formation of the first kind of the Prandtl 
secondary flows. In turbulent flow the large-scale vor- 
ticity is generated by the divergence of the Reynolds 
stresses. This mechanism determines the second kind 
of the Prandtl turbulent secondary flows ■ 

The generation of large-scale vorticity in a homoge- 
neous nonhelical turbulence with an imposed large-scale 
linear velocity shear has been recently studied in 
Let us discuss a mechanism of this phenomenon. The 
equation for the mean vorticity W = VxU read 



9W 



Vx(UxW + F-iyVxW) 



(1) 



where U is the mean fluid velocity, F.; — j {uiUj) is 
the effective force caused by velocity fluctuations, u, and 
V is the kinematic viscosity. The first term, UxW, in 
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Eq. IT]) determines laminar effects of the mean vorticity 
production caused by the sheared motions, while the ef- 
fective force F determines the turbulent effects on the 
mean fluid flow. Let us consider a simple large-scale lin- 
ear velocity shear U*^") = (0, Sx, 0) imposed on the small- 
scale nonhelical turbulence. The equation for the pertur- 
bations of the mean vorticity, W — {Wx{z),Wy{z),{)), 
reads 



dt 
~dt 



SWy + iy^W^ , 



(see [ni), where W" = d^W/dz^ 



(2) 
(3) 

is the turbu- 
lent viscosity, Ig is the maximum scale of turbulent mo- 
tions and the parameter /3o is of the order of 1, and de- 
pends on the scaling exponent of the correlation time 
of the turbulent velocity field (see Sect. II). A solu- 
tion of Eqs. (O and ^ has the form cx exp{'jt + iK^z), 
where the growth rate of the large-scale instability is 
given by 7 = \/]3o S Iq Kz — v^. and Kz is the wave 
number. The maximum growth rate of perturbations of 
the mean vorticity, 7niax = Po {S Iq)'^ / Av^ , is attained at 
Kz — Km = Sl{)/2v^. This corresponds to the ratio 
Wy/Wx = VPqIq Km ~ Stq, where the time tq = Iq/uq 
and mq is the characteristic turbulent velocity in the max- 
imum scale Iq of turbulent motions. Note that in a lam- 
inar fiow this instability does not occur. 

The mechanism of this instability is as follows (see [ll| 
for details). The first term, SWy = (W^^) • V) Ox, in 
Eq. ((2) determines a " skew- induced" generation of per- 
turbations of the mean vorticity Wx by stretching of the 
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equilibrium mean vorticity W(") = (0,0,5), where U 
are the perturbations of the mean velocity. In partic- 
ular, the mean vorticity W^e^ is generated from WyBy 
by equilibrium shear motions with the mean vorticity 
W("), whereby W^e^ oc (W^'^ ■V)U^e^ cx Wyey x W^"). 
Here e^;, and are the unit vectors along x, y and 
z axes, respectively. On the other hand, the first term, 
—PqSIqW^, in Eq. ^ determines a "Reynolds stress- 
induced" generation of perturbations of the mean vor- 
ticity Wy by the Reynolds stresses. In particular, this 
term is determined by (VxF)^. This implies that the 
component of the mean vorticity WyBy is generated by 
an effective anisotropic viscous term cx — ^ {WxSx • 
V)U^''\x)ey (X —llSW^'ey. This instability is caused 
by a combined effect of the sheared motions ("skew- 
induced" generation) and the "Reynolds stress- induced" 
generation of perturbations of the mean vorticity. 

The mechanism for this large-scale instability in a 
sheared nonhelical hornogeneous turbulence is different 
from that discussed in [12, [13, [l3| , where the generation 
of large-scale vorticity in the helical turbulence occurs 
due to hydrodynamic alpha effect. The latter effect is 
associated with the hydrodynamic helicity of turbulent 
flow. In a nonhelical homogeneous turbulence this effect 
does not occur. 

The large-scale instability in a nonhelical homogeneous 
turbulence has been studied in [TT] only for a simple case 
of unbounded turbulence with an imposed linear velocity 
shear and when the perturbations of the mean vorticity 
depend on one spatial variable z. In this study the theo- 
retical approach proposed in is further developed and 
applied for comprehensive investigation of the large-scale 
instability for different situations with nonuniform shear, 
inhomogeneous turbulence and a more general form of 
the perturbations of the mean vorticity W(r) that de- 
pends on three spatial variables. 

In the present study we consider three types of the 
background large-scale flows, i.e., the Couette flow (lin- 
ear velocity shear) and Poiseuille flow (quadratic velocity 
shear) in a small-scale homogeneous turbulence, and the 
"log-linear" velocity shear in an inhomogeneous turbu- 
lence. We have derived new mean-field equations for per- 
turbations of large-scale velocity which depend on three 
spatial coordinates in a small-scale sheared turbulence, 
for a nonuniform background large-scale velocity shear 
and for an arbitrary scaling of the correlation time T(fc) 
of the turbulent velocity field. 

The stability of the laminar Couette and Poiseuille 
flows in a problem of transition to turbulence has been 
studied in a number of publications (see, e.g., [l^, [H, 
[13, [H, [S [lOl , and references therein). It is known that 
laminar plane Couette flow and antisymmetric mode of 
laminar plane Poiseuille flow are stable with respect to 
small perturbations for any Reynolds numbers. A sym- 
metric mode of laminar plane Poiseuille flow is stable 
when the Reynolds number is less than 5772 In 
laminar flows the ToUmien-Schlichting waves can be ex- 
cited. The molecular viscosity plays a destabilizing role 



in laminar flows which promotes the excitation of the 
ToUmien-Schlichting waves (see, e.g., [IBl)- These waves 
are growing solutions of the Orr-Sommerfeld equation. 

In the present study we have found a turbulent ana- 
logue of the ToUmien-Schlichting waves. These waves 
are excited by a small-scale sheared turbulence, i.e., by a 
combined effect of the turbulent Reynolds stress-induced 
generation of perturbations of the mean vorticity and the 
background sheared motions. The energy of these waves 
is supplied by the small-scale sheared turbulence. We 
demonstrate that the off-diagonal terms in the turbulent 
viscosity tensor play a crucial role in the excitation of 
the turbulent ToUmien-Schlichting waves. These waves 
can be excited even in a plane Couette flow imposed on a 
small-scale turbulence when perturbations of velocity de- 
pend on three spatial coordinates. When perturbations 
of large-scale velocity depend on one or two spatial co- 
ordinates the turbulent ToUmien-Schlichting waves can 
not be excited in a sheared turbulence. In the present 
study we show that the large-scale Couette and Poiseuille 
flows imposed on a small-scale turbulence can be unstable 
with respect to small perturbations. The critical effective 
Reynolds number (based on turbulent viscosity) required 
for the excitation of this large-scale instability, is of the 
order of 200. 

This paper is organized as follows. In Sect. II the gov- 
erning equations are formulated. In Sect. Ill we consider 
a homogeneous turbulence with a large-scale linear ve- 
locity shear (Couette flow), while in Sect. IV we study a 
homogeneous turbulence with a large-scale quadratic ve- 
locity shear (Poiseuille flow). In Sect. V we investigate 
formation of large-scale vortical structures in an inhomo- 
geneous turbulence with an imposed nonuniform velocity 
shear. Finally, we draw conclusions in Sec. VI. 



II. GOVERNING EQUATIONS 

The equation for the mean velocity U in incompressible 
flow reads 

+ U • [/, = - ^ + V, {u,uj) + i^AU, , (4) 

where U is the mean velocity, P is the mean pressure 
and 1/ is the kinematic viscosity. The effect of turbulence 
on the mean flow is determined by the Reynolds stresses 
{uiUj), where u are the fluid velocity fluctuations. 

We consider a turbulent flow with an imposed mean 
velocity shear ViU*^*-', where U'"-*. In order to study 
a stability of this equilibrium we consider perturbations 
U of the mean velocity, i.e., the total mean velocity is 
U = U^*' -I- U. Thus, the linearized equation for the 
small perturbations of the mean velocity is given by 

(^| + u(^).v)[7, + (U.V)t/i^) = -^ + F. 

+ yAU^, (5) 
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where Fi = —Vj /ij(U) is the effective force, /y = {uiUj) 
and P are the perturbations of the fluid pressure. Equa- 
tion ^ is derived by subtracting Eq. (|4]) written for the 
equiUbrium velocity U^"' from Eq. ^ for the mean veloc- 
ity U. We consider a simple large-scale velocity shear, so 
that U'-''-* is directed along y direction and is non-uniform 
in X direction, i.e., U^'*) = {0,U^'\x),0). 

In order to obtain a closed system of equations, an 
equation for the effective force Fi = — Vj/y (U) has been 
derived in where 



+C3 Gij] , 



(6) 



{dU)ij = {S7iUj + VjUi)/2 and Iq is the maximum scale 
of turbulent motions. The tensors M^-, iY^j, Hij and 
Gij, in the expression for the Reynolds stresses ([6]) are 
given by: 



M,, = {dU^%^{dU)„rj+{dU^'^)j„,{dU)r 



Hi- 

Gi; 



WJe. 



n ['-■ntrn 



(at/)™, -I- , 



-■njiri 



Eijk is the fully antisymmetric Levi-Civita tensor, 
(at/('*))jj = (y^U^/^ +VjU^"^)/2 and the parameters Cfc 
in Eq. ([6|) are given below. 

The effective force Fi depends on the correlation time 
of the turbulent velocity field T(fc), where k is the wave 
number. In the present study we derive a more general 
form of the effective force Fi for an arbitrary scaling of 
the correlation time r(fc) — C tq {k/ko)^'^ of the tur- 
bulent velocity field, where ko = I/Iq. To this end we 
use Eq. (20) derived in [ill . The value of the coefficient 
C = {q — 1 + ^)/{q — 1) corresponds to the standard 
form of the turbulent viscosity in the isotropic turbu- 
lence, i.e., Vj. = J T{k) [(u^) E{k)] dk = ro(u^)/3. Here 
E{k) = (q—l) k^^ (fc/fco)^'' is the energy spectrum of tur- 
bulence. For the Kolmogorov's type background turbu- 
lence (i.e., for the turbulence with a constant energy flux 
over the spectrum) , the exponent ^ — q — 1 and the co- 
efficient C = 2. This case has been studied in For a 
turbulence with a scale-independent correlation time, the 
exponent /x = and the coefficient C = 1. The parame- 
ters Cfc entering in the Reynolds stresses ([5]) are given by 
Ci = 2C^ (^2 _ 11^ + 28)/315, C2 = -C2 {7fi+ l)/90 
and C3 = -C2 (^ + 3)/90. 

For the derivation of the effective force Fi we use a 
procedure outlined below (see [ll| for details). Using the 
equation for fluctuations of velocity written in a Fourier 
space, we derive equation for the two-point second-order 
correlation function of the velocity fluctuations (uiUj). 
We introduce a background turbulence with zero gradi- 
ents of the mean fluid velocity. This background turbu- 
lence is determined by a stirring force that is independent 
of gradients of the mean velocity. In this study we use 



a model of isotropic, homogeneous and nonhelical back- 
ground turbulence. Then we subtract the equation for 
the two-point second-order correlation function of the ve- 
locity fluctuations {uiUj)^^'> written for the background 
turbulence from the equation for {ui Uj). This yields the 
equation for the deviations from the background turbu- 
lence. 

The obtained second-moment equation include the 
first-order spatial differential operators N applied to the 
third-order moments M'^^^^\ A problem arises how to 
close the equation, i.e., how to express the third-order 
terms J\fM^^^^'> through the lower moments Aff^^^^ (see, 
e-g-, HH, m, nil). To this end we use a spectral r 
approximation which postulates that the deviations of 
the third-moment terms from the contributions to these 
terms afforded by the background turbulence are ex- 
pressed through the similar deviations of the second mo- 
ments (see, e.g., [n!, [H], [13, [H, i!] ) . A justification of 
the T approximation for different situations has been per- 
formed in numerical simulations and analytical studies in 

[13, HI m, 113, [Mi. 

We assume that the characteristic time of variation of 
the second moment of velocity fluctuations is substan- 
tially larger than the correlation time for all turbulence 
scales. This allows us to obtain a steady state solution of 
the second moment equation for the deviations from the 
background turbulence. Integration in k space allows us 
to determine the Reynolds stresses in the form of Eq. ([6]). 
Note that this form of the Reynolds stresses in a turbu- 
lent flow with a mean velocity shear can be obtained even 
by simple symmetry reasoning (see [ll| for details). 

In the next Sections we use Eq. ([5]) with the derived 
effective force (see Eq. ([S])) for a study of the dynamics 
of perturbations of the mean velocity. We show that 
under certain conditions the large-scale instability can 
be excited which causes formation of large-scale vortical 
structures. 



III. LINEAR VELOCITY SHEAR (COUETTE 
FLOW) IN HOMOGENEOUS TURBULENCE 

We consider a homogeneous turbulence with a mean 
linear velocity shear, U'^'*^ — (0, Sx, 0). This velocity field 
is a steady state solution of the Navier-Stokes equation. 
Let us first study the case when the velocity perturba- 
tions U(i, X, z) are independent of y. The equations for 
the components Ux and Uy of the velocity perturbations 
read 



d_ 

di 



A 



dt 



A 



Uy = -SAUx , 



(7) 
(8) 



and the component Uz is determined by the continu- 
ity equation V-U = 0, where (3o ^ Ci + C2 — C3 = 
(2//2 - 43 + 63) /315. In order to derive Eqs. ^ and 
^ we calculate Vx (VxU) using Eq. ([5]), that allows us 
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to exclude the pressure term from this equation. We also 
use Eq. ^ for the Reynolds stresses in the sheared tur- 
bulence. For simplicity, in Eq. ([8]) we neglect the small 
terms ^ O[{lo/Ls)'^], where Ls is the characteristic scale 
of the velocity shear. 

We seek for a solution of Eqs. ^ and fS]) in the form 



Ux,y = exp(7t) [Ax.y cos{Kx x) 
X cos{Kz z + (p) , 



Bx,y cos\i{Kz x) 



(9) 



where the coefhcients A^^y, Bx.y, the angle 



and the 

growth rate 7 of the instability are determined by the 
boundary conditions. We choose the symmetric solution 
(relative the point a; = 0), because the maximum growth 
rate of the symmetric mode is higher than that of anti- 
symmetric mode (see below). Perturbations of the mean 
velocity grow in time due to the large-scale instability 
with the growth rate 



(10) 



The maximum growth rate of perturbations of the mean 
velocity, 



7n 



(11) 



is attained at — Km = \/]3o S lo/2v^. 

In order to determine the threshold required for the 
excitation of the large-scale instability, we consider the 
solution of Eqs. ([7]) and ^ with the following bound- 
ary conditions for a layer of the thickness Ls in the x 
direction: at a; = ztLs/2 the functions U = and 
{Ux,y) — 0. This yields the threshold value of the 
wave number K", determined by the equation 

tan(i^^'-Ls/2) = - tanh(i^^'X5/2) . (12) 

The condition 7,nax > implies that Km > K". There- 
fore, the large-scale instability is excited when the value 
of the shear S exceeds the critical value Scr that is given 

by 



<5'cr Tq 



4.7 



Ls ' 



(13) 



where i^T"' = 2tt/Ls- Note that the value of K" for 
the the symmetric mode is smaller than that for anti- 
symmetric mode. This is the reason why the maximum 
growth rate of the symmetric mode is larger than that of 
antisymmetric mode. 

Note that the parameter (3o depends on the scaling ex- 
ponent /i of the correlation time of the turbulent velocity 
field, T(fc) c>c k^^. In particular, for the Kolmogorov scal- 
ing, T(fc) oc /c"^/"^, we arrive at (3o — 0.45. This case has 
been considered in [ij]. The necessary condition for the 
large-scale instability (/So > 0) reads 2 /i^ — 43 -I- 63 > 0, 
i.e., the instability is excited when < /i < 1.58 and 
fi > 19.9. Note that the condition /i > 19.9 is not realis- 
tic. In the case of a turbulence with a scale-independent 



correlation time, the exponent /i = and the parameter 
/3o = 0.2. 

For small hydrodynamic Reynolds numbers, the scal- 
ing of the correlation time T(/e) ~ l/{i'k'^), i.e., fi — 2, 
and the parameter /3o < 0. This implies that the in- 
stability of the perturbations of the mean vorticity does 
not occur for small Reynolds numbers in agreement with 
the recent results obtained in ^32i] whereby an instability 
of the perturbations of the mean vorticity in a random 
flow with large-scale velocity shear has not been found 
using the second order correlation approximation and as- 
sumption that the correlation time T{k) ~ l/{vk'^). This 
approximation is valid only for small Reynolds numbers 
(see discussion in [s^l)- 

Let us consider now a more general case when the 
velocity U depends on three spatial coordinates, i.e., 
U = U(i, x, z). The equations for the components Ux 
and Uy of the velocity perturbations read 



_9 

at 



dt 



+ V,, ~ Aj ^Ux ^llS ^ [/3o ^HUy 

+ (/3i-/?2)V,V^J7,] , (14) 

+ C/(^) V, -u^l^Uy^llSl^ [132 (A - V^) Ux 

+ (/?i - /3o) VxVy Uy] + S{2Vl -~A)Ux, (15) 



and the component Uz is determined by the continu- 
ity equation V-U = 0. Here A^ = A - V^, /3i = 
2C1-C2 = C2 (8Ai2-39^-H231)/630and/?2 = C1+C3 = 
(4^2 _ 51^^ _|_ 9l)/630. In order to derive Eqs. ^ 

and wc calculate Vx(VxU) using Eq. (O, that 
allows us to exclude the pressure term from this equa- 
tion. For the derivation of Eqs. and ([TS]) we also 
use Eq. ([6]) for the Reynolds stresses in the sheared tur- 
bulence. Equations and can be reduced to the 
Orr-Sommerfeld equation if we replace vt by v and set 
/?„ = (see, e.g., [H, [H, [13 j and references therein). 

We seek for a solution of Eqs. p4)l - (|15p in the form 
(X ^{x) exp{jt + iLLjt + i Kh ■ r), where K-h is the wave 
number that is perpendicular to the a;-axis. After the 
substitution of this solution into Eqs. p^ - p^ we obtain 
the system of the ordinary differential equations which is 
solved numerically. We consider the solution of Eqs. p4|) - 
P^ with the following boundary conditions for a layer 
of the thickness Ls in the x direction: at a; = ztLs/2 
the functions U = and V^, {Ux,y) = 0. These boundary 
conditions with a linear velocity shear corresponds to the 
Couette flow. 

In this Section we show that in a small-scale turbulence 
the large-scale Couette flow can be unstable under cer- 
tain conditions. The range of parameters {Ls/Lh', f) for 
which the large-scale instability occurs is shown in Fig. [U 
where Lr = 2t: /Kh, Kh = {K^ + K'^)^/'^ and if is the 
angle between the wave vector K/f and the direction of 
the mean sheared velocity U^*-'. In Figs. [2111] we show 
the growth rate of the large-scale instability 7 Tq and the 
frequencies of the generated modes utq versus Ls/Lh- 
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FIG. 1: Range of parameters [Ls / Lh', for which the large- 
scale instability occurs for Couette background flow and for 
different values of the large-scale shear: S to — 0.2 (dashed 
line), Sto = 0.4 (solid line). Here Ls/lo = 30. 



The growth rates of the large-scale instability increase 
with the increase of the angle ip, while the frequencies of 
the generated modes decrease with the angle (p so that 
uj{(p 90°) 0. The growth rate of the large-scale 
instability reaches the maximum value at if = 90°. In 
addition, the range of angles ip for which the large-scale 
instability occurs, is small and located in the vicinity 
of if = 90° (see Fig. [1]). Therefore, Ky < and since 
Lz ^ Ls, the size of the structures in the direction of U^^-' 
is much larger than the sizes of the structures along x 
and z directions. This implies that the large-scale struc- 
tures formed due to this instability are stretched along 
the mean sheared velocity U'^"). 

The curves in Figs. [2111] have a point whereby the 
first derivative of the growth rate of the large-scale in- 
stability with respect to the wave number Kh has a sin- 
gularity. At this point there is a bifurcation which is 
illustrated in Fig. [31 In particular, the growth rates and 
the frequencies for the first and the second modes which 
have the highest growth rates are shown in Fig. [5^ and[3jD. 
When the size of perturbations Lh < i*, the frequen- 
cies of the first and the second modes are different, but 
the growth rates are the same. Therefore, at the point 
Lh = i*, there is a generation of two different modes 
with the same growth rate. On the other hand, when the 
size of perturbations Lh > i*, the growth rates of the 
first and the second modes are different, but the frequen- 
cies are the same. 

The maximum growth rate of perturbations of the 
mean velocity, 7max7 is attained at Kh = Km, and the 
value Km increases with the increase of the angle ip be- 
tween the wave vector K.h and the direction of the mean 
sheared velocity U^*^ . The increase of shear S promotes 
the large-scale instability, i.e., it cause the increase of the 
range for the instability (see Fig. [1]) and the maximum 
growth rate (see Figs. [21 and [J]). The characteristic spa- 
tial scale Lm = ^n/Km and the time scale tinst ~ "f, 
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FIG. 2: The growth rate (a) of the large-scale instability 
7 To and frequencies ro of the generated modes (b) ver- 
sus Ls/Lh for Couette background flow and for different 
angles if.if = 85° (dashed-dotted line), = 87° (dashed 
line), ip = 90° (solid line). Here S to — 0.4, Ls/lo ~ 30 and 
io{(p = 90°) = 0. 



for the instability are much larger than the character- 
istic turbulent scales. This justifies separation of scales 
which is required for the validity of the mean-field theory 
applied in the present study. The spatial profiles of the 
ratios of vorticity components Wy/Wx and Wz/Wx for 
perturbations in Couette background flow are shown in 
Fig. [3 The function Wy/Wx is symmetric relative to the 
center of the flow at x = 0, while the function Wz/Wx 
is antisymmetric. Since the function Wx — > at the 
boundaries of the flow, the ratios of vorticity components 
Wy/Wx and Wz/Wx tend to ^ ±oo at the boundaries. 

The numerical results for the case ip = 90° shown in 
Figs. [21 [H and [5l coincide with the analytical predictions 
based on Eqs. (|9|)- (fT3ll . For instance, the threshold value 
of the shear at Ls/lo = 30 is SctTq « 0.157 in agree- 
ment with Eq. (|13p . The ratio of vorticity components 



1 

max 



Wy/Wx ~ 0.3 at X = for modes with the maximum 
growth rate of the large-scale instability. This is in agree- 
ment with this ratio of Wy/Wx obtained using Eq. 
The maximum growth rates of perturbations of the mean 
velocity are in agreement with Eqs. ([TT]) and (I12|) . When 
we switch off the turbulence, the large-scale instability 
does not excited, etc. 

The growing modes with a nonzero frequency discussed 



6 




0.0 0.2 0.4 0.6 0.8 1.0 L^/Lj 




FIG. 3: The growth rate 7ro (a) and the frequency uitq (b) 
versus Ls / Lh of the first (sohd hne) and the second (dashed 
hne) modes which have the highest growth rates for Couette 
background flow. Here the angle <^ = 87°, S to = 0.4 and 
Ls/k = 30. 




FIG. 4: The growth rate (a) of the large-scale instability 
7 TO and frequencies uj to of the generated modes (b) ver- 
sus Ls/Lh for Couette background flow and for different 
angles i^: = 87° (dashed-dotted line), <^ = 88° (dashed 
line), = 90° (solid line). Here Sto ^ 0.2, Ls/lo = 30 and 
Lo{ip = 90°) = 0. 



in this Section can be regarded as the turbulent ana- 
logue of the ToUmien-Schhchting waves. In laminar flows 
the ToUmien-Schlichting waves are growing solutions of 
the Orr-Sommerfeld equation and the molecular viscos- 
ity promotes the excitation of the ToUmien-Schlichting 
waves (see, e.g., [3). On the other hand, the turbulent 
Tollmien-Schlichting waves are excited by a small-scale 
sheared turbulence, i.e., by a combined effect of the tur- 
bulent Reynolds stress-induced generation of perturba- 
tions of the mean vorticity and the background sheared 
motions. 



IV. QUADRATIC VELOCITY SHEAR 
(POISEUILLE FLOW) IN HOMOGENEOUS 
TURBULENCE 



Now we consider a homogeneous turbulence with an 
imposed large-scale quadratic velocity shear, U^^^ = 
S<,x{\ — x/Ls)ey. The equations for the components 
Ux and Uy of the velocity perturbations read 



{PoAh 



(/3i - /?2) V, ?7J + S' Vy , (16) 



A]Uy = S{2S/l - A) a 



-2S'V.,Ux+llSA iPl- Pa) VccVyUy 



lis' 



+132 (A - V^) 

~VyUx) + A [{2^2 + Pi) y.u. 

+ iP2-PQ)^yUy 



2piVxVy{V^Uy 



(17) 



and the component Uz is determined by the continu- 
ity equation V-U = 0, where S{x) = V x U'^'^^ and 
S' = VxS. In order to derive Eqs. IHl) and we 
calculate Vx(VxU) using Eq. ([5]). We seek for a solu- 
tion of Eqs. and P?|) in the form oc ^(x) exp(7i -|- 
iujt + i K/f • r), where Kjj is the wave number that is per- 
pendicular to the a;-axis. After the substitution of this 
solution into Eqs. and ([T7)) we obtain the system 
of the ordinary differential equations which is solved nu- 
merically. We consider the solution of Eqs. P^ - (fT7|) with 
the following boundary conditions for a layer of the thick- 
ness Ls in the x direction: at a; = ± Ls/2 the functions 
U = and [Ux.y) = 0. These boundary conditions 
with a quadratic large-scale velocity shear corresponds to 
the Poiseuille flow. We show below that in a small-scale 
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FIG. 5:_ The spatial profiles of the ratios of vorticity compo- 
nents Wy/Wx (a) and Wz/Wx (b) for modes with the max- 
imum growth rates of the large-scale instability in Couette 
background flow and for different angles ip: ip — 85° (dashed- 
dotted line), ip = 87° (dashed line), ip = 90° (solid line). Here 
Sto = 0.4 and Ls/lo = 30. 
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FIG. 7: The growth rate (a) of the large-scale instability 
7 To and frequencies ui tq of the generated modes (b) ver- 
sus Ls/Lh for Poiseuille background flow and for different 
angles ip: ip — 84° (dashed-dotted line), ip = 87° (dashed 
line), (p = 90° (solid line). Here S* to = 0.6, Ls/lo = 30 and 
uj{(p = 90°) = 0. 




0.5- 

O.OI ■ , ■ : ' : ' , ' 

80 82 84 86 88 (p 



FIG. 6: Range of parameters {Ls/Lh; v) for which the large- 
scale instability for Poiseuille background flow occurs, and for 
different values of the large-scale shear: S, to = 0.5 (dashed 
Une) and S* to = 0.6 (solid line). Here Ls/lo ~ 30 and S* — 
S{x = 0). 



turbulence the large-scale Poiseuille flow can be unstable 
with respect to small perturbations. 

The range of parameters {Ls/Lh] ^p) for which the 
large-scale instability in the Poiseuille background flow 
occurs is shown in Fig. [S] for different values of the large- 
scale shear, where S"* = S{x = 0). The growth rates of 



this instability and the frequencies of the generated tur- 
bulent ToUmien-Schlichting waves are shown in Figs. [7| 
and [5] The spatial profiles of the ratios of vorticity com- 
ponents Wy/Wx and Wz/Wx in Poiseuille background 
flow for modes with the maximum growth rates of the 
large-scale instability are shown in Fig. [5] The general 
behaviour of the large-scale instability in the Poiseuille 
background flow is similar to that for the Couette back- 
ground flow. In particular, the growth rates of the large- 
scale instability increase with the increase of the angle 
tp between the wave vector K/f and the direction of the 
mean sheared velocity U^"^ , reaching the maximum value 
a,t (p = 90°. The frequencies lutq of the generated turbu- 
lent ToUmien-Schlichting waves by the large-scale insta- 
bility decrease with the increase of the angle ip and ui ^ 
a.t (f ^ 90°. The values Km at which the growth rates 
of the large-scale instability reach the maximum values 
increase with the increase of the angle ip. The range for 
the large-scale instability and the growth rates of pertur- 
bations in the Poiseuille background flow increases with 
the increase of shear. This implies that increase of shear 
promotes the large-scale instability. 

For the Poiseuifle flow the large-scale instability can 
be excited for smaller angles (p than that for the Cou- 
ette background flow. On the other hand, the thresholds 
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FIG. 8: The growth rate (a) of the large-scale instability 
7 To and frequencies uj tq of the generated modes (b) ver- 
sus Ls/Lh for Poiseuille background flow and for different 
angles ip: tp = 87° (dashed-dotted line), ip — 88° (dashed 
line), ip = 90° (solid line). Here 5** ro = 0.5, Ls/lo = 30 and 
uj{ip = 90°) = 0. 
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FIG. 9: The spatial profiles of the ratios of vorticity compo- 
nents Wy/Wx (a) and Wz/Wx (b) for modes with the max- 
imum growth rates of the large-scale instability in Poiseuille 
background flow and for different angles ip: ip = 84° (dashed- 
dotted line), tp = 87° (dashed line), (p = 90° (solid line). Here 
5. TO = 0.6 and Ls/lo = 30. 



for the instability in the value of shear and in the value 
of Ls/Lh for Poiseuille background flow are larger than 
that for the Couette background flow. A difference be- 
tween the Couette and Poiseuille background flows can 
be also seen in Figs. [5] and [5] for the spatial profiles of 
the ratios of vorticity components Wy/Wx and Wz/Wx- 
This difference is caused by the different geometries in 
these flows. In particular, the first spatial derivatives of 
the fiow velocity in the Poiseuille background flow are 
antisymmetric relative to the center of the fiow at a; = 0, 
while they are symmetric (constant) in the Couette back- 
ground flow. This is the reason of that the spatial profile 
of Wy/Wx is symmetric relative to x = in the Couette 
background fiow, and it is antisymmetric in the Poiseuille 
flow. 



V. NONUNIFORM VELOCITY SHEAR IN 
INHOMOGENEOUS TURBULENCE 

In this Section we consider a more complicated form 
of nonuniform velocity shear in an inhomogeneous turbu- 
lence. For simplicity we consider the case when the small 
perturbations of the mean velocity U are independent of 
y. The equations for the components Ux and Uy of the 



velocity perturbations in an inhomogeneous turbulence 
with a nonuniform shear read 
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VlUy 



2{\7l„]\/lUx , 



d_ 

di 



.A 



Ux^f3o llSA-Vl{llS) 

(18) 

-S + I3i Vx{llS)Vx 

V,i^^) V,C/j, , (19) 



Uy = 



Ux 



and the component Uz is determined by the continu- 
ity equation V-U = 0, where S{x) = V^, C/^*^ Equa- 
tion (fT9|) is the y component of Eq. ^ with Vj,P = 0, 
while Eq. (jTS]) is the x component of Vx(VxU) de- 
termined from Eq. ([5]). We consider the solution of 
Eqs. (fT8|) and (fT9|) with the following boundary condi- 
tions for a layer of the thickness Ls in the x direction: 
at a; = ± Ls/2 the functions U = and Wx {Ux.y) = 0. 

We consider a " log- linear" velocity proflle for the back- 
ground large-scale flow in an inhomogeneous turbulence. 
In particular, we use the following relationship for the 
velocity shear S{x) = u1/vrp{x) and the eddy viscosity 
v^{x) = u^^Iq{x), where lo{x) = Kr]{x)Ls is the turbu- 
lence length scale, k is the von Karman constant, Ui, is the 
friction velocity, r]{x) is the dimensionless function that 




0.0 0.2 0.4 0.6 0.8 1.0 v;(x) 

FIG. 10: The spatial profile of the normahzed turbulent vis- 
cosity v^{x) for different values of the parameter a: a = 6 
(solid line), a — 10 (dashed line), a = 20 (dotted line), a = 50 
(dashed-dotted line). 
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FIG. 11: The mean velocity profile U'-"' (x) /Um^^ for diflterent 
values of the parameter a: a — 6 (solid line), a = 10 (dashed 
line), a = 20 (dotted line), q = 50 (dashed-dotted line), where 



characterizes the spatial profile of the background ve- 
locity shear and inhomogeneity of small-scale turbulence 
(see below). These relationships are usually used for the 
logarithmic boundary layer profiles (see, e.g., [2^). The 
spatial profile rj^x) for < x < Ls/2 is chosen in the 
form 

r]{x) — ai [l ~ exp(— flo i)] + a2 x + x'^ + a4 , 

(20) 

where x = x/ Ls — 1/2, the coefficients Ofe are determined 



FIG. 12: The growth rate of the large-scale instability versus 
imax/iz in inhomogcneous turbulence with nonuniform veloc- 
ity shear for different values of the parameter a: a = & (solid 
line), a = 10 (dashed line), q = 20 (dotted line) and a = 50 
(dashed-dotted line). Here imax = lo{x ^ 0.5 Lz). 



by the following conditions: at x the functions rj — 1, 
V^Tj = 0, V^77 = 0, Vlrj = 0, and at x = ~Ls/2 the 
derivative VxV — ct/Ls- Here a is a free parameter that 
characterizes the inhomogeneities of small-scale turbu- 
lence. The spatial profile of the normalized turbulent 
viscosity h'*{x) = v^{x)/{KUi, Ls) = ri{x) is shown in 
Fig. [To] for different values of the parameter a. The func- 
tion v* (x) is chosen to be symmetric relative the point 
a; = 0. The minimum possible value of the parameter 
a is a = 6. We have chosen the velocity shear profile 
[/^''-'(a;) so that the logarithmic velocity profile near the 
boundaries can be matched with the linear shear velocity 
for the central part of the background flow. Such kind of 
flow is typical for the atmospheric boundary layer. Fig- 
ure [11] shows the mean velocity profile U'^'^\x)/U-nia,x for 
different values of the parameter a, where C/max — u^,/ n- 

We seek for a solution of Eqs. (|18p and in the form 
cx '^{x) exp(7i -I- iKzz). After the substitution of this 
solution into Eqs. (US]) and we obtain the system of 
the ordinary differential equations which is solved numer- 
ically. The growth rate ^tq of the large-scale instability 
versus ^max/^z is shown in Fig. ll2i where — 27: / is 
the size of perturbations in z direction and /max = k is is 
the maximum value of the turbulent length scale Iq when 
77 ^ 1 {x ^ \). The range of parameters (/max/^z; 
a) for which the large-scale instability occurs is shown 
in Fig. [T5b . The vertical dashed line in Fig. [T3] indi- 
cates that the minimum possible value of the parame- 
ter a is amin = 6. Figure [T5b demonstrates that the 
increase of the parameter a causes the increase of the 
maximum growth rate of the large-scale instability. The 
growth rate of the large-scale instability for the inhomo- 
geneous turbulence with a large-scale nonuniform shear 
is much larger than that for the Couette and Poiseuille 
background flows. 

The spatial profiles of the ratios of vorticity compo- 
nents Wy/Wx and Wz/Wx for modes with the maximum 
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FIG. 13: (a). The range of parameters (imax/i/z; a) for which 
the large-scale instability in inhomogeneous turbulence with 
nonuniform velocity shear occurs, (b). The maximum growth 
rate 7max tq of the large-scale instability versus the parameter 
a. Here Zmax = lo{x 0.5 L^). 
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FIG. 14: The spatial profiles of the ratios of vorticity com- 
ponents Wy/Wx (a) and Wz/Wx (b) for modes with the 
maximum growth rates of the large-scale instability in in- 
homogeneous turbulence with nonuniform velocity shear for 
different values of the parameter a: a = 6 (solid line), 
Of = 20 (dashed line) and a = 50 (dashed-dotted line). Here 



lo(x — + 0.5 1/z 



growth rates of the large-scale instability are shown in 
Fig. [131 These profiles are different from that for the 
Couette and Poiseuille background flows. The compo- 
nents Wy and Wz of perturbations of the mean vorticity 
in the central part of the flow are usually much smaller 
than the component Wx- Inspection of Figs. [T^ and [T^ 
shows that the parameter /max/iz < 0.17. The charac- 
teristic time scale for the instability is much larger than 
the characteristic turbulent time. This justifles separa- 
tion of scales which is required for the validity of the 
mean-field theory used here. 

Note that in the interval —Ls/2 < x < the obtained 
results discussed in this Section imply a stability theory 
for the turbulent boundary layer. Our study shows that 
the turbulent boundary layer can be unstable under cer- 
tain conditions. 



VI. DISCUSSION 

In this study the theoretical approach proposed in [llj 
is further developed and applied to investigate the large- 
scale instability in a nonhelical turbulence with a nonuni- 
form shear and a more general form of the perturba- 
tions of the mean vorticity. In particular, we consider 
three types of the background large-scale sheared flows 



imposed on small-scale turbulence: Couette flow (lin- 
ear velocity shear) and Poiseuille flow (quadratic veloc- 
ity shear) in a small-scale homogeneous turbulence, and 
a more complicated nonuniform velocity shear with the 
logarithmic velocity profile near the boundaries matched 
with the linear shear velocity for the central part of the 
background fiow. This nonuniform velocity shear is im- 
posed on an inhomogeneous turbulence. The latter fiow 
is typical for the atmospheric boundary layer. 

We show that the large-scale Couette and Poiseuille 
flows imposed on a small-scale turbulence are unstable 
with respect to small perturbations due to the excita- 
tion of the large-scale instability. This instability causes 
generation of large-scale vorticity and formation of large- 
scale vortical structures. The size of the formed vorti- 
cal structures in the direction of the background velocity 
shear is much larger than the sizes of the structures in 
the directions perpendicular to the velocity shear. There- 
fore, the large-scale structures formed during this insta- 
bility are stretched along the mean sheared velocity. In- 
crease of shear promotes the large-scale instability. The 
thresholds for the excitation of the large-scale instability 
in the value of shear and the aspect ratio of structures 
for Poiseuille background flow are larger than that for the 
Couette background flow. The growth rate of the large- 
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scale instability for the inhomogeneous turbulence with 
the "log- linear" velocity shear is much larger than that 
for the Couette and Poiseuille background flows. The 
characteristic spatial and time scales for the instability 
are much larger than the characteristic turbulent scales. 
This justifies separation of scales which is required for the 
validity of the mean-field theory applied in the present 
study. 

The large-scale instability results in excitation of the 
turbulent Tollmien-Schlichting waves. The mechanism 
for the excitation of these waves is different from that 
for the Tollmien-Schlichting waves in laminar flows. In 
particular, the molecular viscosity plays a crucial role in 
the excitation of the Tollmien-Schlichting waves in lam- 
inar flows. Contrary, the turbulent Tollmien-Schlichting 
waves are excited by a combined effect of the turbulent 
Reynolds stress-induced generation of perturbations of 
the mean vorticity and the background sheared motions. 
The energy of these waves is supplied by the small-scale 
sheared turbulence, and the off-diagonal terms in the tur- 
bulent viscosity tensor play a crucial role in the excitation 
of the turbulent Tollmien-Schlichting waves. 

Note that this study is principally different from the 
problems of transition to turbulence whereby the stabil- 
ity of the laminar Couette and Poiseuille flows are in- 
vestigated (see, e.g., [H [11 [H [H, [H, US], and refer- 
ences therein). Here we do not analyze a transition to 
turbulence. We study the large-scale instability caused 
by an effect of the small-scale anisotropic turbulence 
on the mean flow. This anisotropic turbulence is pro- 
duced by an interaction of equilibrium large-scale Cou- 
ette or Poiseuille flows with a small-scale isotropic back- 
ground turbulence produced by, e.g., a steering force. 
The anisotropic velocity fluctuations are generated by 
tangling of the mean- velocity gradients with the velocity 
fluctuations of the background turbulence [III HI] . 

The "tangling" mechanism is an universal phenomenon 
that was introduced in [H, [sH] for a passive scalar and 
in [H 1131 for a passive vector (magnetic field). The 
Reynolds stresses in a turbulent flow with a mean veloc- 
ity shear is another example of tangling anisotropic fluc- 



tuations [38[. For instance, these velocity fluctuations are 
anisotropic in the presence of shear and have a steeper 
spectrum oc k^'^l'^ than, e.g., a Kolmogorov background 
turbulence (see, e.g., [2a [H, [S^, |40, |4l| ) . The anisotropic 
velocity fluctuations determine the effective force and the 
Reynolds stresses in Eq. This is the reason for the 
new terms c>c /3„ Zq appearing in Eqs. (fTi|) - P^ . 

The obtained results in this study may be of rele- 
vance in different turbulent astrophysical, geophysical 
and industrial flows. Turbulence with a large-scale ve- 
locity shear is a universal feature in astrophysics and 
geophysics. In particular, the analyzed effects may be 
important, e.g., in accretion disks, extragalactic clusters, 
merged protostellar and protogalactic clouds. Sheared 
motions between interacting clouds can cause an excita- 
tion of the large-scale instability which results in gener- 
ation of the mean vorticity and formation of large-scale 
vortical structures (see, e.g., [i3j[S[ll|)- E)ust particles 
can be trapped by the vortical structures to enhance ag- 
glomeration of material and formation of particle clusters 

14l[4i[i3,[ii,[ii. 

The suggested mechanism can be used in the analy- 
sis of the flows associated with Prandtl's turbulent sec- 
ondary flows (see, e.g., p?, ^To'|). However, in this study 
we have investigated only simple physical mechanisms to 
describe an initial (linear) stage of the formation of vor- 
tical structures. The simple models considered in this 
study can only mimic the flows associated with turbu- 
lent secondary flows. Clearly, the comprehensive numer- 
ical simulations of the nonlinear problem are required for 
quantitative description of the turbulent secondary flows. 
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